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ABSTRACT 

We devise a Monte-Carlo based, optimized filter match method to extract the 
thermal Sunyaev-Zel'dovich (SZ) signature of a catalog of 116 low-redshift X-ray 
clusters from the first year data release of the Wilkinson Microwave Anisotropy 
Probe (WMAP). We detect an over-all amplitude for the SZ signal at the ~ 8cr 
level, yielding a combined constraint of /gas^ = 0.08 ± O.Ol(ran) ± 0.01{sys) 
on the gas mass fraction of the Intra-Cluster Medium. We also compile X-ray 
estimated gas fractions from the literature for our sample, and find that they are 
consistent with the SZ estimates at the 2a level, while both show an increasing 
trend with X-ray temperature. Nevertheless, our SZ estimated gas fraction is 
30 — 40% smaller than the concordance ACDM cosmic average. We also express 
our observations in terms of the SZ flux-temperature relation, and compare it 
with other observations, as well as numerical studies. 

Based on its spectral and spatial signature, we can also extract the microwave 
point source signal of the clusters at the 3a level, which puts the average mi- 
crowave luminosity (at ~ 41 GHz) of bright cluster members [M^ < ~21) 
at (2.4 ± 0.8) X 10^^/i~^erg/s/Hz . Furthermore, we can constrain the aver- 
age dark matter halo concentration parameter to Cvir = 3.4lQ g, for clusters with 

> 5 keV. 

Our work serves as an example for how correlation of SZ surveys with cluster 
surveys in other frequencies can significantly increase our physical understanding 
of the intra-cluster medium. 
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1. Introduction 

Clusters of galaxies are the largest relaxed concentrations of mass in the universe. They 
are interesting for cosmology as they probe the evolution of the large scale structure of the 
universe (e.g. Eke, Cole, & Frenk 1996; Viana & Liddle 1999; Haiman, Mohr, & Holder 2001; 
Hu & Haiman 2003) and they are interesting on their own, as we can resolve, probe, and 
study their inner structure in different frequencies, ranging from microwave to X-rays, and 
also through weak and strong gravitational lensing of background galaxies (e.g., Carlstrom, 
Holder, & Reese 2002; Nichol 2004; Markevitch 2003; Hennawi & Spcrgcl 2004; Massey et al. 
2004; Sand, Treu, Smith, & Ellis 2004). What adds to this simplicity is that, at least for the 
massive clusters, almost all of the baryonic matter sits in the diffuse ionized Intra-Cluster 
Medium (ICM), which can be studied theoretically and observationally with relatively simple 
physics, and give us a census of cosmic baryonic budget (see e.g.. White, Navarro, Evrard, 
& Prenk 1993; Evrard 1997; Mohr et al. 1999). 

In this paper, we focus on the microwave signatures of galaxy clusters, the thermal 
Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 1972), caused by the scattering of 
Cosmic Microwave Background (CMB) photons by hot gas in the diffuse ICM, and yielding 
characteristic spatial and spectral imprints on the CMB sky. 

The thermal SZ effect has changed from the subject of theoretical studies to that of 
intense observational endeavor within the past decade, as various experiments have and are 
being designed to study this effect (e.g., APEX, ACT, AMI, Planck, SZA, SPT; see Bond 
2002, for an overview). The main reason behind this wide attention is the potential for 
using SZ detected clusters as standard candles to probe the cosmological evolution up to 
large redshifts (e.g., Haiman, Mohr, & Holder 2001; Verde, Haiman, & Spergel 2002; Lin & 
Mohr 2003; Majumdar & Mohr 2003). Compared to the SZ detection method, the X-ray 
detected clusters, which have been primarily used for this purpose until now (e.g., Henry 
2000; Vikhlinin et al. 2003; Henry 2004), become much harder to detect at large redshifts, 
and are also believed to be more affected by complex astrophysics associated with galaxy 
formation, cooling or feedback within clusters (Carlstrom, Holder, & Reese 2002). 

Moreover, unlike X-ray observations which only sample regions of high gas density, 
thermal SZ observations probe the distribution of thermal energy in the cluster, and thus 
provide independent information about the over-all thermal history (Is there an entropy 
floor?; e.g., Voit, Bryan, Balogh, & Bower 2002; Ponman, Sanderson, & Finoguenov 2003), 
and baryonic budget of the cluster (Are there missing baryons?; e.g., Cen & Ostriker 1999). 

Although various scaling relations of X-ray properties of clusters have been extensively 
studied in the literature, mainly due to the scarcity, incoherence, or low sensitivity of SZ 
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observations of clusters, there have been only a few statistical analyses of SZ scaling prop- 
erties in the literature (Cooray 1999; McCarthy, Holder, Babul, & Balogh 2003; Benson et 
al. 2004). Given the upcoming influx of SZ selected cluster catalogs, a good understanding 
of these scaling relations, and in particular, the SZ flux-Mass relation (sec §3.1), which is 
of special significance for cosmological interpretations of these catalogs (e.g. Majumdar & 
Mohr 2003), is still lacking. 

While the first year data release of the Wilkinson Microwave Anisotropy Probe (WMAP; 
Bennett ct al. 2003a) has constrained our cosmology with an unprecedented accuracy, due 
to its low resolution and low frequency coverage, the SZ effect cannot be directly observed 
in the WMAP CMB maps (Huffenberger, Seljak, & Makarov 2004). One possible avenue is 
cross-correlating CMB anisotropies with a tracer of the density (which traces clusters and 
thus SZ signal) in the late universe (Peiris & Spergel 2000; Zhang & Pen 2001). In fact, 
different groups have reported a signature of anti-correlation (which is what one expects 
from thermal SZ at WMAP frequencies) at small angles between WMAP maps and different 
galaxy or cluster catalogs, at a few sigma level (Bennett et al. 2003b; Fosalba, Gaztanaga, 
& Castander 2003; Fosalba & Gaztanaga 2004; Myers et al. 2004; Afshordi, Loh, & Strauss 
2004). While thermal SZ is the clear interpretation of this signal, relating such observations 
to interesting cluster properties can be confused by the physics of non-linear clustering or 
galaxy bias. 

Hernandcz-Montcagudo & Rubiho-Martm (2004) and Hernandcz-Montcagudo, Genova- 
Santos, & Atrio-Barandcla (2004) use an alternative method, where they construct SZ tem- 
plates based on given cluster or galaxy catalogs, and then calculate the over-all amplitude 
of WMAP signal temperature decrement associated with that template. While the method 
yields significant SZ detections (2-5(t), the physical interpretation is comphcated by the 
non-trivial procedure that they use to construct these templates. 

In this paper, we follow the second line by devising an optimized filter match method 
based on an analytic model of ICM which is motivated by both numerical simulations and 
observations. We then apply the method to a sample of X-ray clusters, and construct tem- 
plates of both SZ and potential point source contamination based on the X-ray temperatures 
of each cluster. Combining these templates with the WMAP maps yields constraints on the 
physical properties of our ICM model, namely the ICM gas mass fraction, and the dark 
matter halo concentration parameter. 

Almost all the SZ observations up to date use an isothermal (3 model, obtained from 

X-ray observations, to describe the cluster SZ profile. However, it has been demonstrated 
that, as X-rays and SZ cover different scales inside the cluster, such extrapolation can lead 
to errors as big as a factor of 2 in the interpretation of SZ observations (Schmidt, Allen, & 
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Fabian 2004). Instead, Schmidt, Allen, & Fabian (2004) suggest using a physically motivated 
NFW profile (see §3) to model the X-ray and SZ observations simultaneously, which in 
their case, leads to consistent estimates of Hubble constant for three different clusters. We 
choose to follow their approach in choosing a physically motivated ICM model, rather than 
a mathematically convenient one. 

In Appendix A, we introduce a semi-analytic NFW-based model for the ICM gas profile. 
We then start in §2 by describing the WMAP CMB temperature maps and our compiled X- 
ray cluster catalog. §3 derives the theoretical SZ/Point Source templates based on our ICM 
model, while §4 describes oiu statistical analysis methodology. In §5 we describe the results 
of our statistical analysis, listing the constraints on gas fraction, concentration parameter, 
and point source contamination of our clusters. It the end, §6 discusses the validity of various 
assumptions that we made through the treatment, and §7 highlights the major results and 
concludes the paper. 

Throughout the paper, we assume a ACDM flat cosmology with = 0.3, and Hq — 
lOO/i km/s/Mpc. While no assumption for h is made in our analysis of the SZ signal, we 
adopt the value oi h — 0.7 to compare the X-ray gas fractions with our SZ signal. 

2. Data 

2.1. WMAP foreground cleaned CMB maps 

We use the flrst year data release of the observed CMB sky by WMAP for our analysis 
(Bennett et al. 2003a). The WMAP experiment observes the microwave sky in 5 frequency 
bands ranging from 23 to 94 GHz. The detector resolution increases monotonically from 0.88 
degree for the lowest frequency band to 0.22 degree for the highest frequency. Due to their low 
resolution and large Galactic contamination, the two bands with the lowest frequencies, K(23 
GHz) and Ka(33 GHz), are mainly used for Galactic foreground subtraction and Galactic 
mask construction (Bennett et al. 2003b), while the three higher frequency bands, which 
have the highest resolution and lowest foreground contamination, Q(41 GHz), V(61 GHz), 
and W(94 GHz), are used for CMB anisotropy spectrum analysis (Hinshaw et al. 2003). 
Bennett et al. (2003b) combine the frequency dependence of 5 WMAP bands with the 
known distribution of different Galactic components that trace the dominant foregrounds 
(i.e. synchrotron, free-free, and dust emission) to obtain the foreground contamination in 
each band. This foreground map is then used to clean the Q, V and W bands for the 
angular power spectrum analysis. Similarly, we use the cleaned temperature maps of these 
three bands for our SZ analysis. We also use the same sky mask that they use, the Kp2 
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mask which masks out 15% of the sky, in order to avoid any remaining Galactic foreground. 
However, we stop short of masking out the 208 identified WMAP point sources, as many 
of them happen to be close to our clusters. For example, there are 29 WMAP identified 
microwave sources within 3.6 degrees of 66 of our clusters. Instead, we decide to model the 
point source contamination based on its frequency dependence (§3.2). 

The WMAP temperature maps and mask are available in the HEALPix format of spher- 
ical coordinate system (Gorski, Hivon, & Wandelt 1999), which is an equal-area pixeliza- 
tion of the sphere. The resolution of the first year data is Ngide = 512 = 2''. implying 
12 X A^gj^^g = 3, 145, 728 independent data points (in lieu of masks) and ~ 0.1° sized pixels, 
for each sky map. 

2.2. Cluster Catalog and X-ray Data 

Our objective is to study the SZ signal in a large sample of galaxy clusters. To this end, 
we have assembled our sample from several existing X-ray cluster samples (David et al. 1993; 
Mohr et al. 1999; Jones & Forman 1999; Finoguenov, Reiprich, & Bohringcr 2001; Reiprich & 
Bohringer 2002; Ikebe et al. 2002; Sanderson et al. 2003), as X-ray observations may provide 
reliable cluster mass estimates, and avoid false detections due to chance projections. The 
selection criteria require that the clusters (1) must have measured X-ray emission weighted 
mean temperature (Tx), (2) are reasonably away from the Galactic plane (Galactic latitude 
|6| > 10°), and (3) are at least 3 degrees away from the Galactic foreground Kp2 mask (§2.1). 

The redshift information is obtained from NED and/or SIMBAD, and the above cata- 
logs. The cluster temperature is taken from the literature cited above, primarily the study 
of Ikebe et al. (2002). We have adopted the Tx obtained when the central cool core region 
is excluded, and identify the peak of the X-ray emission (either from the cluster catalogs or 
from archival ROSAT images) as the cluster center. Our final cluster catalog contains 117 
nearby clusters, (0 < 2; < 0.18), whose temperature ranges from 0.7 to 11 keV. 

Our requirement that clusters have measured Tx is to provide reliable mass estimates. 
Given Tx, the observed mass-temperature relation (Finoguenov, Reiprich, & Bohringer 2001, 
hereafter FRBOl) 

M500 = (1.781°:?? X 10'=^/i-'Mo) rx(keV)^-5«-o.o? (1) 

can be used to obtain M500 = (47r/3)500pc'"5oo' mass enclosed by rsoo, within which the 
mean overdensity is 500 times the critical density of the universe pc- The wide range of 
cluster temperature in our sample implies that our clusters span two orders of magnitude in 
mass. 
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In order to facilitate comparisons of our SZ-derived gas fraction with the X-ray mea- 
surements, we compile the gas mass (within r5oo) for most of our clusters from the literature 
(Mohr et al. 1999; Jones & Forman 1999), supplemented by the data based on the study of 
Sanderson et al. (2003). The gas masses provided by Jones & Forman (1999) arc measured 
at a fixed metric radius of 0.5h^^ Mpc; we convert it to the nominal radius of rsoo by the 
measured /3-model profile, and then to the virial radius r2oo (= ''"vir; see §3), using the analytic 
model of §3. 

Fig. (1) shows the distribution of redshifts and virial radii for our clusters (using the 
analytic model of §3 for dark matter concentration Cvir = 5). The solid lines show the 
resolution limits of the three WMAP bands, as well as the physical radius of the 1 degree 
circle at the cluster redshift. We see that most of our clusters are in fact resolved in all the 
WMAP bands. 



In Appendix A, based on the assumption of hydrostatic equilibrium and NFW dark 
matter profile (Navarro, Frenk, & White 1997), we develop an analytic model for the gas 
and temperature distribution in the Intra-Cluster Medium (ICM). In this model, assuming a 
given NFW concentration parameter, Cyir, all the properties of the cluster/ICM are quantified 
in terms one parameter, which, can be taken to be e.g., the cluster virial mass, Myir, or its 
X-ray temperature (see Equations A2 and All for definitions). In particular, Mvir can be 
expressed in terms of (the observed) Tx (Equation All)^ within the model for a given value 
of Cvir- Now, let us estimate the dominant microwave signals of a galaxy cluster, based on 
our simple model. 



The contribution of the thermal SZ effect to the CMB temperature anisotropy (see 
Carlstrom, Holder, & Reese 2002, for a review), at the frequency u, is proportional to the 
integral of electron pressure along the line of sight 



-'^Notc that we have assumed the observed X-ray temperature, Tx, to be the emission weighted temperature 
in our model. We address the error introduced due to this assumption in §6.2. 



3. Modeling the Intra-Cluster Medium 



3.1. SZ profile 




(2) 
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where ar is the Thomson scattering cross-section, and nie is the electron mass. The SZ flux, 
deflned as the integral of STsz over the solid angle, is then given by 

STsz -6^=1 STsziiy; n) d'h = _ ^^^cmb^(W^cmb) f p^^^^ . ^-2^^^ ^3^_ 
Jsn "^eC Jsn 

Here, dA is the angular diameter distance, and x spans over the cone extended by the solid 
angle Sfl. Now, assuming local thermal equilibrium 

the total SZ flux of a cluster is 

(1 + X) ar T,,e F{hi^/TcMJ /gas M^^, T^, 



2 nip rrieC^ d\ 



-(1.42 X 10-^ mK) (0.1 deg)^ f ^WT^cmb) /y ^ j ^^^^^j^^y) f 



which can be combined with our analytic model (equations A8, All, and A12) to obtain 

Stotir^) = -(2.16 X 10-^ mK) (0.1 deg)^ 1 TTTTTFTTTTl^ I ^'/'(keV) [B{c^,,) f,^ h] 

L E{z) [HodA{z)/c\ } 

- -(1.63 X 10- mJy) | ■ ,2, ^^11, , w l4 ^'"'^^^^^ ^^^^^^^ 

Lsmh {x/2)E{z)[HodA{z)/c\ ) 

where 

^ = W^CMB. (7) 

^ \ Jg xMx ) ■ V/ P/^ xUx ) ' 

£;(^) = i/(^)/i/o = [^^m(l + ^)' + ^^a] , (9) 
HodA{z)/c^{l + zr'£^^, (10) 

and functions /(x) and g'(a:) are deflned in Appendix A (Equations A7 and AlO). For the 
relevant range of 3 < Cvir < 8, which is consistent with various measurements of cluster dark 
matter profile (see Lin, Mohr, & Stanford 2004 for a brief review), i?(cvir) is a decreasing 
function of Cvir which varies from 2.1 to 1.4. Note that all the factors in equation (6), with 
the exception of the last one, arc fixed by observations. Therefore, [-B(cvii.)/gas/i] is the 
combination of model parameters which will be fixed by our SZ fiux observations. 
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We then use a Monte-Carlo method to reproduce the expected SZ flux of a given cluster. 
In this method, we equally distribute the total SZ flux of equation (6) among Nmc = 4000 
random points, whose 3D density follow the ICM pressure, Pg{r) oc f{r/r^-^^ ■ g{r/r^ij.), 
around the center of a given cluster. While the method would be equivalent to exact projec- 
tion in the limit Nmc ~^ oo, the poisson error introduced due to a finite value of Nmc will 
be negligible comparing to the WMAP detector noise. The projected distribution of points 
should be then smeared by the WMAP beam window to get the expected distribution of the 
SZ flux. The expected SZ signal of the cluster in pixel i is proportional to rij, the number of 
points that will fall into that pixel: 

Si{iy) = 5Tsz,i ■ 5^^pixei = (^] X Stotiiy). (11) 



3.2. Point Source Contamination 

The frequency dependence of WMAP small angle anisotropics have been interpreted 
as a random distribution of point sources with a flat spectrum (i.e. Antenna temperature 
scahng as Bennett et al. 2003b). The majority of individually identifled WMAP point 
sources are also consistent with a flat spectrum. Since the SZ signal has a small frequency 
dependence at the range of WMAP frequencies (41 GHz < i/ < 94 GHz), we can use 
this frequency dependence to distinguish the Point Source (PS) contamination from the SZ 
signal. To do this, we assume a microwave point source with a flat (constant) luminosity per 
unit frequency, Lps, for each cluster galaxy, and that the galaxies follow the dark matter 
distribution (equation Al) inside each cluster. The total microwave flux associated with the 
point sources is then given by 

where A^vir is the number of galaxies, above a certain magnitude limit, within the virial radius, 
and (11 is the luminosity distance. For our analysis, we use the Lin, Mohr, & Stanford (2004) 
result for 2MASS near infrared X-band selected galaxies: 

\ 0.85±0.04 



iV.,(M.< -21) = 37 ±3 ^ ^^^^.3-.^ J . (13) 

Thus, Lps is defined as the total point source luminosity per unit frequency associated 
with the cluster, divided by the number of galaxies brighter than the near infrared i^-band 
magnitude of —21, within the virial radius of the cluster. 
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4. Statistical Analysis Methodology 

For a low resolution CMB experiment such as WMAP, the main sources of uncertainty 
in the SZ signal are the primary CMB anisotropics, as well as the detector noise. Since both 
of these signals are well described by gaussian statistics, we can write down the which 
describes the likehhood of observing a given model of the cluster SZ+PS profile (see §3.1): 

i,j;a,b 

where a and b run over WMAP frequency bands (i.e. Q, V, or W), and i and j run over the 
WMAP pixels. Here, Tj^ and Si{'Ua) are the observed temperature and expected SZ+PS flux 
in pixel i and band a, while Cia,jb is the covariance matrix of pixel temperatures: 

Cia,jl, = nl5i^5a, + Y.{^^^ |W^pixel(^)r W^beam(^; o) VFbeam(^; 6) QP,(cOS%). (15) 

Here, riia is the pixel detector noise, VFpixei and Wbcam are the HEALPix pixel and WMAP 
beam transfer functions (Page et al. 2003), C^'s and P^'s are the primary CMB multipoles 
and Legendre polynomials respectively, and 9ij is the angular seperation between the pixels 
i & j. We use CMBfast code (Seljak & Zaldarriaga 1996) in order to generate the expected 
values of Q's for the WMAP concordance ACDM cosmology (Bennett et al. 2003a). 

Because WMAP detector noise only varies on large angular scales, riia can be assumed 
to be almost constant if we limit the analyses to the neighborhood of a cluster, yielding 




Kab{(-) = W^beam(^; a)W^,,r^{t, h) C, + n^(5„;,(50pixel, (17) 

where we used 

<^^^pixel (^^) |l^pi.el(^)rA(cOS%) = 5,,. (18) 

Now, it is easy to check that, in the small angle limit, we have 

CraJ6 ^ <^^^pixel E |W^pixel(^)r ^aV(^) ^^(cOS%), (19) 

We can again use the Monte-Carlo method, described at the end of §3.1, to evaluate 

C~-'^S = C~-'^WbeamSo, where So is the raw projected SZ profile. To do so, instead of 
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smearing So by the detector beam windoWjWbeam, we can smear So by C ^Wbeam, which is 
given by 

(20) 

Since the SZ signal is dominant at small angles, and at the same time we want to avoid 
the non-trivial impact of the CMB masks on the covariance matrix inversion, we cut off Diajb 
if the separation of pixels i and j is larger than 6d — 3°. We do not expect this to impact 
our analysis significantly, as the CMB fluctuations are dominated by smaller angles. As 
Ti — Si/5Q,pixei still minimizes x^, this truncation cannot introduce systematic errors in our 
SZ or PS signal estimates. However, it may cause an underestimate of the covariance errors. 
In §4.1, we introduce a Monte-Carlo error-estimation method to alleviate this concern. 

The x^, given in equation (14), is quadratic in Ai = B{c)fgash and A2 = Lps, and can 
be re-written, up to a constant, as 

X' = F«/^A„A;3-2A-A„, (21) 

where 

F"' = Sn^L E Sn^a)C.y,S^{u,), (22) 

i,j-.a,b 

and = ^ Sn^a)Cr;jbTjb. (23) 

i,j;a,b 

Note that Sii^Ua) = Ao.S'f (z/^) is the sum of the SZ -|- PS flux contributions per pixel (derived 
in §3.1 and §3.2) for all the clusters in the sample. 

After evaluating the coefficients F"'^ and A" via the Monte-Carlo method described 
above, the x^ in Eq. (21) can be minimized analytically to obtain the best fit values for the 
gas fraction and point source luminosity. After this minimization, the resulting x^ can be 
used to constrain the value of the concentration parameter Cyir- 



4.1. Error Estimates 

While the covariance matrix obtained from the x^ in equation (14) gives a natural way 
to estimate the errors, our Monte-Carlo based approximation of the covariance matrix, as 
well as its truncation beyond On = 3°, may reduce the accuracy of our error estimates. 
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Another source of error, which is not included in the covariance method, is the uncertainty 
in observed X-ray temperatures. 

In order to obtain more accurate error estimates, we use our primary CMB power 
spectrum (from CMBfast), combined with the WMAP noise and beam properties to generate 
99 Monte-Carlo Realizations of WMAP CMB maps in its three highest frequency bands (Q, 
V, & W). Neglecting the contamination of cluster signals by background points sources, 
these maps can then be used to estimate the error covariance matrix for our fgash and Lps 
estimators, within an accuracy of V2/99 14%. 

To include the impact of Tx errors in our Monte-Carlo error estimates, we assume 
an asymmetric log- normal probability distribution for the true temperature, P(Tx), which 
is centered at the observed value, T°^*, and its extent on each side is given by the the 
upper /lower error of the observed temperature, 6T^/6Tl, i.e. 



Therefore, in each Monte-Carlo realization, the temperature of each cluster is also randomly 
drawn from the above distribution, which is then used to construct the SZ/PS template for 
that cluster (§3.1). 



In this section, we use the framework developed in §4 to combine the WMAP tempera- 
ture maps with our cluster catalog. It turns out that about 30% of our clusters are within 
3°, and about 8% within 1° of another cluster in our sample, implying possible correlations 
between the signals extracted from each cluster. However, in order to simplify the analysis 
and interpretation of our data, we ignore such possible correlations, and thus assume that 
the values of /gas and Lps, obtained for each cluster is almost independent of the values for 
the rest of the sample. As we argue below, there is no evidence that this approximation may 
have biased our error estimates of global averages significantly. 

One of our clusters (A426; Perseus chistcr) shows an 18cr {Lps = (1.69 ± 0.09) x 
10^^ /i~^erg/s/Hz/galaxy) signature for frequency dependent PS signal. It turns out that 
the 5th brightest microwave source detected by the WMAP team (WMAP#94; NGC 1275) 
happens to be the brightest galaxy of the cluster. As this point source overwhelms the SZ 
signal, we omit A426 from our analysis, which leaves us with a sample of 116 X-ray clusters. 



where 




(24) 



5. 



Results 
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Table 1: The global best fit values to the gas fraction and point source parameters for different 
temperature cuts, and assuming Cyir — 3 (top) or Cyir = 5 (bottom), compared with X-ray 
estimates (see the text) . Ax^ shows the relative significance of the best fit with respect to 
the no-cluster (null) model. 



Cyir 3 


# 


/gas/i (X-ray) 


/gas/i (SZ) 


Lp5(1028/i-2erg/s/Hz/gal) 




all clusters 


116 


0.0608 ±0.0004 


0.077 ±0.011 


0.19 ±0.08 


-53.1 


> 3 keV 


78 


0.0743 ±0.0005 


0.073 ±0.012 


0.07 ±0.16 


-47.8 


Tx > 5 keV 


38 


0.0954 ±0.0007 


0.086 ±0.015 


0.07 ±0.44 


-55.6 


Tx > 8 keV 


8 


0.1151 ±0.0016 


0.083 ±0.020 


0.06 ±0.49 


-27.2 





# 


/gas/i (X-ray) 


/gas/i (SZ) 


Lp5(1028/i-2erg/s/Hz/gal) 




all clusters 


116 


0.0662 ± 0.0004 


0.084 ±0.011 


0.24 ±0.08 


-57.3 


Tx > 3 keV 


78 


0.0810 ±0.0005 


0.084 ±0.013 


0.21 ±0.21 


-51.2 


Tx > 5 keV 


38 


0.1040 ±0.0008 


0.101 ±0.017 


0.40 ±0.53 


-60.6 


Tx > 8 keV 


8 


0.1254 ±0.0018 


0.094 ± 0.021 


0.33 ±0.59 


-29.8 



5.1. Global ICM gas fraction and Point Source Luminosity 

The most straightforward application of the statistical framework introduced in §4 is to 
obtain a global best fit for the gas fraction and galaxy microwave luminosity Lpg for a 
given value of concentration parameter c^i^. Table 1 shows the results of our global fits for 
nominal values of Cyir = 3 and Cyir = 5, within different temperature cuts, which are also 
compared with the estimates from our compiled X-ray observations. Note that the lower 
value of Cyir is probably appropriate for the high end of the cluster masses/temperatures, 
while the higher value may correspond to less massive clusters. To get the X-ray gas fraction, 
the gas mass estimated from X-ray observations (§2.2) is divided by the virial mass expected 
from observed Tx (Eq. All) for each value of Cyir. 

While the overall significance of our model detections are in the range of 5 — 8(T, we 
see that the significance of our SZ detection is ~ Scr for the whole sample, and there is a 
signature of point source contaminations at ~ 3(7 level, although we should note that there 
is a significant correlation (70 — 80%) between the SZ and PS signals. 

While the SZ signal is mainly due to massive/hot clusters, most of the PS signal comes 
from the low mass/temperature clusters (compare 1st and 2nd rows in each section of Table 
1). Therefore, for the PS signal, the higher concentration value of Cyir = 5 might be closer 
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to reality, putting the average microwave luminosity of cluster members at 

(Lp5)(41 GHz; Mk < -21) = (2.4 ± 0.8) x lO^^/i-^erg/s/Hz/galaxy, (25) 

Surprisingly, this number is very close to the diffuse WMAP Q-band luminosity of Milky 
Way and Andromeda galaxy (Afshordi, Loh, & Strauss 2004), i.e. ~ 2 x 10^'''erg/s/Hz. 

Therefore, assuming that a significant fraction of cluster members have a diffuse emission 
similar to Milky Way, our observation indicates that, on average, nuclear (AGN) activity 
cannot overwhelm the diffuse microwave emission from cluster galaxies. Nevertheless, models 
of microwave emission from faint (radiatively inefficient) accretion flows cannot be ruled out 
(see §6.3). 

As an independent way of testing the accuracy of our error estimates, we can evaluate 
the for the residuals of our global fits for the whole sample (first rows in Table 1). For 
Cvir = 3 and Cvir = 5, the residual for our global fits are 252 and 268, respectively, 
which are somewhat larger than (but within 2(7 of) the expected range for 2 x 116 degrees 
of freedom, i.e. 230 ± 22. While this may indicate ~ 7% underestimate of errors, it may 
at least be partly due to the Tx dependence of /gas, which we discuss in the next section. 
Since correlation of errors among close clusters may decrease this value, while systematic 
underestimate of errors tends to increase the residual x^, we conclude that, unless these two 
effects accidentally cancel each other, we do not see any significant evidence (i.e. > 10%) 
for either of these systematics. Repeating the exercise for the hotter sub-samples of Table 1 
yields a similar conclusion. 

Finally, we note that the X-ray and SZ values for fg^sh are always consistent at the 2(7 
level (see further discussion below). 

5.2. Dependence on the Cluster Temperature 

Let us study the dependence of our inferred cluster properties on the cluster X-ray 
temperature, which can also be treated as a proxy for cluster mass (Eq. All). Since 
the errors for individual cluster properties are large, we average them within 2 keV bins. 
The binned properties are shown in Figs. (2) & (3), and listed in Table 2. Similar to the 
previous section, we have also listed estimated gas fractions based on our compilation of 
X-ray observations. 

Fig. (2) compares our SZ and X-ray estimated gas fractions. The solid circles show our 
SZ observations, while the triangles are the X-ray estimates. 

We notice that, similar to the global averages (Table 1), our SZ signals are more or 
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Table 2: The global best fit values to the gas fraction and point source parameters for 
different temperature bins, assuming Cvir = 3 (top) or Cyir = 5 (bottom), compared with X- 
ray estimates (see the text) . Ax^ shows the relative significance of the best fit with respect 
to the no-cluster (null) model. Note that the last bin only contains one cluster, i.e. A2319. 



Tx(keV) 


Tx(kcV) 


# 


/gas/i (X-ray) 


fg^h (SZ) 


Lp<j(1028/;,-2erg/s/Hz) 


Ax' 


0-2 


1.1 


20 


0.0139 ±0.0056 


0.256 ±0.167 


0.22 ±0.10 


-5.0 


2-4 


3.3 


44 


0.0485 ± 0.0004 


0.030 ±0.036 


0.13 ±0.20 


-0.7 


4-6 


4.7 


28 


0.0740 ±0.0010 


0.030 ± 0.030 


-0.05 ± 0.54 


-1.9 


6-8 


6.5 


16 


0.0914 ±0.0009 


0.100 ±0.031 


-0.10 ±1.31 


-27.1 


8-10 


8.6 


7 


0.1126 ±0.0018 


0.076 ± 0.025 


-0.05 ±0.51 


-16.5 


10-12 


11.0 


1 


0.1271 ± 0.0039 


0.133 ±0.053 


2.89 ±2.91 


-11.8 


T,(keV) 


Tx(keV) 


# 


fgash (X-ray) 


/gas/i (SZ) 


Lp5(1028/i-2erg/s/Hz) 


Ax' 


0-2 


1.1 


20 


0.0152 ±0.0061 


0.299 ±0.169 


0.29 ±0.11 


-6.3 


2-4 


3.3 


44 


0.0528 ±0.0005 


0.014 ±0.041 


0.19 ±0.26 


-0.6 


4-6 


4.7 


28 


0.0806 ±0.0011 


0.037 ±0.035 


-0.12 ±0.72 


-2.4 


6-8 


6.5 


16 


0.0996 ±0.0010 


0.123 ±0.036 


0.33 ± 1.62 


-29.9 


8-10 


8.6 


7 


0.1227 ±0.0020 


0.085 ± 0.028 


0.14 ±0.63 


-16.7 


10-12 


11.0 


1 


0.1385 ±0.0043 


0.148 ±0.051 


3.98 ±3.24 


-14.5 



less consistent with the X-ray gas estimates. The for the difference of X-ray and SZ bins 
are 6.7 and 7.6 for Cvir = 3 and c^ij. = 5 respectively, which are consistent with the 1-a 
expectation range of 6 ± 3.5, for 6 random variables. Therefore, we conclude that we see no 
signature of any discrepancy between the SZ and X-ray estimates of the ICM gas fraction. 

Another signature of consistency of our X-ray and SZ data points is the monotonically 
increasing behavior of gas fraction with T^^. 

Indeed, this behavior has been observed in previous X-ray studies (e.g., Mohr et al. 
1999; Sanderson et al. 2003), and has been interpreted as a signature of preheating (Bialek, 
Evrard, & Mohr 2001) or varying star formation efficiency (Bryan 2000). A power law fit to 
our binned data points yields 

/gas/i = (0.069 ± 0.014) ° ' (for Cvir = 3), (26) 



^This trend is also responsible for the fact that a global fit (constant /gas) to the sample with T > 3 keV 
is less significant than a global fit to the smaller sample with T^> 5 keV (see Table 1). 
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and /gas/i = (0.077 ± 0.014) L^^^y ) ° ' (for Cvir = 5), (27) 

where the uncertainties in the normahzation and power are almost un-correlated . 

Fig. (3) shows that, after removing A426 (T^, = 6.4 keV), none of our bins show more 
than 2-a signature for point sources. The fact that the observed amphtude of point sources 
changes sign, and is consistent witli zero implies that any potential systematic bias of the 
SZ signal due to our modeling of the point sources (§3.2) must be negligible. 



5.3. SZ flux- Temperature relation 

Given a perfect CMB experiment, and in the absence of primary anisotropics and fore- 
grounds, in principle, the SZ flux is the only cluster property that can be robustly measured 
from the CMB maps, and does not require any modeling of the ICM, while any measure- 
ment of the gas fraction would inevitably rely on the cluster scaling relations and/or the 
assumption of a relaxed spherical cluster. WMAP is of course far from such a perfect CMB 
experiment. Nevertheless, we still expect the SZ flux measurements to be less sensitive to the 
assumed ICM model (§3), compared to our inferred gas fractions. Therefore, here we also 
provide a SZ flux-Tx scaling relation which should be more appropriate for direct comparison 
with other SZ observations and hydro-simulations. Plugging Eqs. (26) and (27) into Eq. (6) 
yields 

^vir(Jy) dlih-' Mpc) E{z) - -(2.41±0.49) x 10^ L{x) (^^^^) (for Cvir = 3), (28) 

and 



6.6 keV 



^+0.8 
-0.1 



,Svir(Jy) dl(h-^ Mpc) E{z) = -(2.04±0.37) x 10=^ L(x) { —^^7 ) (fo^ ^vir = 5), (29) 



where 



m = , (30) 



a;^(a;coth(a;/2) -4) 
sinh2(a;/2) 

and X = hv /RbT^^^ is the detector frequency in units of the CMB temperature. For the three 
highest frequencies of WMAP, Q(41 GHz), V(61 GHz), and W(94 GHz), L{x) = 3.8,7.6, 
and 13.7 respectively. The flts should hold within 3 keV < Tx < 11 keV, which is the range 
of cluster temperatures which contribute the most to our SZ detection. Notice that the 
difference between the normalizations inferred for two concentrations is comparable to the 
measurement errors. 
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Benson et al. (2004) is the only other group which expresses its SZ observations in terms 
of SZ flux-temperature relation. Our result is consistent with their observations, within the 
relevant temperature range {T^ ~ 9 keV for their sample) . This is despite the higher median 
redshift of their sample {z ~ 0.2 — 0.8), which may indicate no detectable evolution in the 
SZ flux-temperature normalization. Their scaling with temperature, however, is signiflcantly 
shallower than our measurement (oc T^^^'^^^'^), which is in contrast with our scaling (cx T^'^), 
at more than 2a level. This is most hkely due to the difference in the range of temperatures 
that are covered in our analysis. Indeed, the clusters in our three highest temperature bins 
whose temperature coincides with that covered in Benson et al. (2004) show a much shallower 
dependence on temperature (see Fig. 2), which is consistent with their results. 

We should note that we have to use our ICM model of §3 to convert 5'vii. = 5*200 to a flux 
within a much smaller area, 5'25oo, which is reported in Benson et al. (2004). The conversion 
factor is 0.30 and 0.42 for Cyir = 3 and 5 respectively. 

As to comparison with numerical simulations of the ICM, even the most recent studies 
of the SZ effect in galaxy clusters (da Silva, Kay, Liddle, & Thomas 2004; Diaferio et al. 
2004) include only a handful of clusters above — 5 keV. This is despite the fact that 
most observational studies of the SZ effect, including the present work, are dominated by 
clusters with Tx > 5 kcV. Therefore, a direct comparison of our observed SZ fluxes with 
numerical studies is not yet feasible. Instead, we can compare the SZ fluxes for clusters 
around = 5 keV, where the temperature range of observed and simulated clusters overlap. 
Making this comparison, we sec that, for the few simulated clusters with Tx ~ 5 keV in da 
Silva, Kay, Liddle, & Thomas (2004) the SZ fluxes are in complete agreement with our 
observations. However, clusters of Diaferio et al. (2004) are underluminous in SZ by close 
to an order of magnitude. Indeed, Diaferio et al. (2004) also notice a similar discrepancy 
with the SZ observations of Benson et al. (2004). The fact that the results of Diaferio et al. 
(2004) are inconsistent with other simulations and observations, may be indicator of a flaw 
in their analysis. 



5.4. Constraining the Concentration Parameter 

It is clear that the assumption of constant concentration parameter, Cyir, which we have 
adopted up to this point, is an oversimplification. The average value of the concentration 
parameter is known to be a weak function of the cluster mass in CDM simulations (oc M~^'^; 
e.g., NFW, Eke, Navarro, & Steinmetz 2001); even for a given mass, it follows a log-normal 
distribution (Bullock et al. 2001; Dolag et al. 2004), which may also depend on mass (Afshordi 
& Cen 2002). 
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As discussed at the end of §4, we can repeat our Monte-Carlo template making procedure 
for different values of Cvir, which yields quadratic expressions for x^, and thus enables us to 
(after marginalizing over Lps) draw likelihood contours in the f^ash — Cvir plane. From 
Tabic (1), wc sec that most of our SZ signal is due to the clusters with Tx > 5 keV (see 
also Table 2), while Cvir(oc T~^'^^) is expected to stay reasonably constant in this range. 
Therefore, we restrict the analysis to this sample. Fig. (4) shows the result, where we have 
explicitly computed the for integer values of Cvir, and then interpolated it for the values 
in-between. The sohd contours show our 68% and 95% likehhood regions (Ax^ — 2.2 and 
6.2). We see that our data can constrain the concentration parameter of the dark matter 
halos to Cvir = 3.4tQ;g (median ± 68% percent hkelihood). The dotted contours show the 
same likelihoods expected for clusters hotter than 5 keV in the WMAP ACDM concordance 
model, where fgash = 0.hh/Q„i = 0.12 ± 0.01 is the upper limit expected from the WMAP 
concordance cosmology (Spcrgcl et al. 2003), while the range of c^ij. is based on an extension 
of the top-hat model (Afshordi & Cen 2002) for the same cosmology, averaged over the 
masses'^ of our cluster sample, and inversely weighted by the square of temperature errors. 

We see that, while the mean gas fraction is about 30% (~ 2a) smaller than the cos- 
mic upper limit (see the 2nd row in Table 1), our inferred constraint on Cyir is completely 
consistent with the ACDM prediction. 

6. Discussions 

6.1. ICM Gas Fraction 

In §5, we demonstrated that, while we are able to detect the thermal SZ effect at the 
7-8(7 level from the first year data release of WMAP temperature maps, the inferred gas 
fractions are typically smaller than the X-ray estimates, as well as the cosmological upper 
hmit (= Qbh/^m — 0.12±0.01). The SZ (as well as X-ray) observations have been often used, 
in combination with the nucleosynthesis bound on fib, to constrain flm, through replacing 
the upper limit by equality (Myers et al. 1997; Mason, Myers, & Readhead 2001; Grego et al. 
2001; Lancaster et al. 2004). Nevertheless, similar to our finding, such determinations have 
consistently yielded lower values than, the now well-established, upper limit. In fact, cooling 
and galaxy formation do lead to a depletion of the ICM gas. To make the matters worse, 
supernovae feedback can make the gas profile shallower, also leading to smaller baryonic 



^The relation between cluster masses and temperatures (Eq. All) is a function of Cyir itself, but the 
uncertainty introduced in cluster mass estimates as a result, only slightly affects the obtained concentration 
range. 
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fraction within a given radius. After all, clusters may not be such accurate indicators of the 
baryonic census in the universe. 

We note that these processes will affect low mass clusters more strongly than high mass 
ones, and therefore it is natural to expect the massive clusters to be better representative 
of the cosmic baryonic content. Indeed, X-ray-based determinations of the ratio of cluster 
gas to virial mass (for massive clusters) have seemingly been more successful in reproducing 
the cosmic average (after ~ 10% correction for stars; e.g., Lin, Mohr, & Stanford 2003). 
However, we should note that, as the X-ray emissivity is proportional to the square of local 
plasma density, any smooth modeling of the ICM which may be used to infer the gas mass 
from the X-ray map of a cluster, tends to overestimate this value due to the contribution 
of unresolved structure to the X-ray emissivity (i.e. (n^) > (n)^). The hydrodynamical 
simulations can suffer from the same problem, and thus fail to estimate the full magnitude 
of the effect. 

There are several factors which can account for the discrepancy between the gas fractions 
determined here and those estimated from previous X-ray studies. Firstly, since direct 
detection of X-rays from the ICM is rarely possible near the virial radius, a significant degree 
of extrapolation is required to infer gas properties at r2oo- X-ray studies typically assume 
a /5-model form for the gas distribution, with an empirically-motivated index parameter of 
/9 ~ 2/3 (e.g., Jones & Forman 1999), implying /)gas oc at large radii. By contrast, 
our physically-motivated model for the gas density (Eq.A9) yields pgas oc at large radii, 
which produces a lower gas mass within r2oo- Secondly, the dark matter concentration is 
likely to be higher than both values assumed here in less massive halos, as a consequence 
of hierarchical formation. This underestimation of Cvir correspondingly underestimates /gas. 
Thirdly, the effects of non-gravitational heating on the ICM in cooler clusters can act to 
displace the gas beyond the radius where we observe it directly (Mohr et al. 1999; Sanderson 
et al. 2003). Consequently, extrapolating to r2oo/^5oo based on the resulting lower-density 
gas that is observed will lead to an underestimate of the total gas mass. 

The difference in logarithmic slope of Pgas('") large radius between a /9-model with 
P 2/3 and Eq. (A9) is partly due to the effects of non-gravitational physics biasing the 
gas distribution with respect to the dark matter. However, there is some evidence that f3 
may increase at larger radius in some (Vikhlinin, Forman, & Jones 1999), although not all 
(Sanderson et al. 2003) clusters, which would reduce the discrepancy. However, satisfactory 
resolution of this issue will require mapping of the gas distribution out to the virial radius 
in a representative sample of clusters; a task which is complicated by the large angular size 
and low surface brightness of the outer regions of the ICM. 
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6.2. Model Uncertainties 

The relationship between the observed SZ flux and the ICM gas fraction (Eq. 6), rehes 
on the accuracy of the (electron) virial temperature- mass relation (Eq. A12). Although our 
X-ray temperature- mass relation (Eq. All) is consistent with observations, the (emission 
weighted) X-ray temperature only probes the inner parts of the cluster, and there can still 
be significant deviations from our simple picture of the uniform ICM in the cluster outskirts. 

For example, Voit et al. (2003) argue that, compared to a uniform homogeneous accre- 
tion, inhomogeneous accretion will inevitably lead to smaller entropy production. Although 
this may not significantly affect the central part of a cluster, it can significantly change the 
boundary condition (see Eq. A5) behind the accretion shock. This can be also interpreted 
as incomplete virialization which can lead to smaller virial temperatures. 

Let us estimate how much error the model uncertainties may introduce to our mea- 
surements. Neglecting the contribution of radio sources (see §6.3), we can divide the model 
uncertainties into the SZ profile shape, and SZ flux uncertainties. 

In §5, we saw that assuming Cyir = 5 instead of Cvir = 3 may result in ~ 10% difference 
in the inferred gas fraction. With the exception of merging clusters, given the low resolution 
of WMAP maps, this is the level of error that we expect may be introduced due to the 
uncertainty in the proflle shape. 

As to the SZ flux uncertainty in our model (Equation 6), we note that the inferred Tx 
dependence of the SZ flux hinges upon the accuracy of our X-ray temperature-mass relation 
(Equation All). Parameterizing this as 

M,i, = ATl/\ (31) 

in Appendix A, we argue that the systematic uncertainty/error in ^ is ~ 10%, i.e. l^A/A ~ 
0.1. As the total SZ flux is proportional to the estimated virial mass, we also get 

- ^ ~ 0.1. (32) 

Moreover, the virial radius of the cluster, which is crucial in our matched filter method, is 
modified according to 

rvir oc Mis = AV3rV2, (33) 

yielding 

Arvir ACvir 1 /A>1 



0.03. (34) 
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Stotd\ 



The SZ template is then given by 

S{e)^\^^^^jj:{edA/r.u), (35) 
where E(a;) is the normahzed SZ (integrated pressure) profile in our model^: 



!lf{y)9{y)- 



ydy 



E(x) = ^ ^/y[^_ (3g^ 

2/0 f{y)9{y)y^dy 



Therefore, the systematic error in the estimated /gas is given by: 



gas 



, 1 / dlnE 
-I + -{2 + 



3 V dlnx 



^ . _^ . _0.I. (37) 



where we used Equation (34), and the fact that din J] /dlnx ^ —2 in cluster outskirts, where 
most of the SZ signal comes from. Therefore, we see that the expected level of systematic 
error due to the theoretical uncertainty in the SZ profile/flux is ~ 15%, which is comparable 
to our random error. However, the similar systematic error in our constraints on Cyir (~ 3%; 
Equation 34), is significantly smaller than the associated random error. 

Further complication may be introduced by the fact that, as a result of ICM inhomo- 
geneities and incomplete frequency coverage, the observed X-ray temperature, T^, could be 
different from the emission weighted temperature. Tew, derived in Appendix A. Recently, 
Mazzotta et al. (2004) defined a spectral-like temperature, Tgi, which can approximate the 
observed to better than a few percent, and Mazzotta et al. (2004) claim that, in hydro- 
dynamic N-body simulations, T^/ can be lower than T^w by as much as 20-30%. However, 
applying the definition of Tgi to the analytic model of Appendix A (which is clearly less 
structured than both simulations and observations), we find that T,/ is smaller than T^yj by 
only ~ 5%. Therefore, we can ignore the impact of this discrepancy in our analyses. 

Finally, another possibility is the breakdown of Local Thermal Equilibrium (LTE) in 
the cluster outskirts. While the hydrodynamic shocks heat up the ions instantly, the char- 
acteristic time for heating up the electrons (via Coulomb interactions) can be significantly 
longer (Fox & Loeb 1997; Chieze, Alimi, & Teyssier 1998; Takizawa 1999), and thus the 
electron temperature can be lower by as much as 20% in the outer parts of clusters. Since 
the thermal SZ effect is proportional to the electron temperature, the breakdown of LTE 
can be a source of low SZ signals. Interestingly, the effect is expected to be bigger for more 
massive clusters which have longer Coulomb interaction times. 



Note that the gas pressure is proportional to f{x)g{x); see Appendix A. 
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6.3. Radio Source spectrum 

While we assumed an exactly flat spectrum for all our point source contamination, a 
more realistic model would include a random spread in the spectral indices a of the point 
sources. In fact, although the average spectral index of the WMAP identified sources is 
zero (flat; Bennett et al. 2003b), there is a spread of ~ 0.3 in a. of individual sources. 
Moreover, although most of the bright microwave sources have an almost flat spectrum, an 
abundant population of faint unresolved sources may have a different spectral index, and yet 
make a significant contribution to the cluster microwave signal. For example, a VLA survey 
finds that for sources with flux S* > 0.1 Jy, the mean spectral index distribution (within 
8 < z/ < 90 GHz) can be described by a Gaussian whose mean and dispersion are —0.37 and 
0.34, respectively (Holdaway, Owen, & Rupen 1994). At fainter flux limit (~ 20 mJy), the 
CBI experiment finds that the mean index (from 1.4 to 31 GHz) is —0.45, with maximum 
and minimum indexes being 0.5 and —1.32, respectively (Mason et al. 2003). 

In order to test the sensitivity of our SZ signal to the point source spectrum, we repeat 
the analysis for the spectral power indices of —1 and 1. Our over- all SZ signal changes by 
less than 2%, implying the insensitivity of our results to the assumed spectrum. 

The average luminosity in the microwave band of cluster galaxies also does not sen- 
sitively depend on our choice of the spectral index. For example, assuming a = —0.8, as 
most studies in low frequencies adopt (e.g., Cooray et al. 1998), we find that at 41 GHz the 
luminosity changes less than 10% to {Lps) = (2.2 ± 0.8) x 10^'''/i~^erg/s/Hz/galaxy. 

Finally, it is interesting to compare our inferred galaxy luminosity at 41 GHz with 
that observed at lower frequencies. A recent study of a large sample of nearby clusters has 
determined the cluster AGN luminosity function (LF) at 1.4 GHz (Lin & Mohr 2005), which 
is in agreement with the bivariate LF obtained by Ledlow & Owen (1996); it is found that the 
cluster LF is very similar to that of the field, once the difference in the overdensity has been 
taken into account. Integrating the LF from 10^^ to 10^^-^ erg/s/Hz (corresponding to the 
observed luminosity range of AGNs) and multiplying by the cluster volume gives the total 
luminosity L1.4 at 1.4 GHz. For a 5 x 10^"^ M0 cluster, L1.4 = 2.43 x 10^^ erg/s/Hz, which 
can be compared to the total luminosity inferred from our result Ltot = (-^P5)A^Mr(l-4/41)" 
(c.f. Eqn 13). Assuming a = 0, we find L^t = 7.12 x 10^^ erg/s/Hz; with a = —0.8, 
Ltot = 9.72 X 10^° erg/s/Hz. 

This exercise suggests that a typical spectrum of a(1.4, 41) ~ — 1 brings our results into 
good agreement with the 1.4 GHz measurements. However, this does not necessarily imply 
the spectral shape to be similar within the frequency range that is relevant to our analysis 
(41 to 94 GHz). The fact that changing a(41, 94) from to -0.8 affects (Lps) less than 10% 
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suggests our analysis is robust against specific choices of the spectral index. 

An alternative to this picture is drawn in Pierpaoli & Perna (2004), where it is proposed 
that a significant part of the point source microwave emission in sky may come from faint (ra- 
diatively inefficient) accretion flows around supermassive black holes in early-type galaxies. 
The assumed microwave luminosity per early- type galaxy in their model is ~ lO^^erg/s/Hz, 
which is consistent with our results, if less than 50% of cluster galaxies with Mk < — 21 are 
early-type galaxies with microwave luminosities at this level. Using the luminosity function 
of near infrared galaxies, we find that the number density of early-type galaxies (Kochanek 
et al. 2001) with Mk < —21 (~ 7.6 x 10~^/i^ Mpc~^) is indeed close to the number den- 
sity of early-types assumed in Pierpaoli & Perna (2004), which is 5.8 x 10^'^/i'^ Mpc^^. A 
nearby optical study of galaxy population in clusters suggests that the early-type fraction 
in clusters is ~ 20%, and does not show strong variation with cluster-centric distance (Goto 
et al. 2003). Based on this finding, we compute the average early-type fraction for galaxies 
with Mk < —21 from the type-specified luminosity functions of Kochanek et al. (2001), 
and assume that it roughly stays the same within cluster virial radii. The obtained fraction 
is ~ 34% (< 50%), which suggests that our observed microwave point source luminosity is 
consistent with the characteristics of the faint accretion fiows proposed in Pierpaoli & Perna 
(2004). 



7. Conclusions and Future Prospects 

In this work, using a semi-analytic model of the Intra-Cluster Medium, we devised a 
Monte-Carlo based optimal filter match method to extract the thermal SZ signal of identified 
X-ray clusters with measured X-ray temperatures. We apply the method to a catalog of 
116 low redshift X-ray clusters, compiled from the literature, and detect the SZ signal, at 
~ 8(7 level (random error), while we estimate a comparable systematic error due to model 
uncertainties. We also see a 3(7 signature for point source contamination, which we model 
based on the assumed spectral characteristics and spatial distribution of point sources. It 
turns out that the average luminosity of bright cluster members {Mk < —21) is comparable 
to that of Milky Way and Andromeda. While our observed SZ signal constrains the gas 
fraction of the Intra-Cluster Medium to 60 — 70% of the cosmic average, it is completely 
consistent with the gas fractions based on our compiled X-ray gas mass estimates. 

Based on our results, we also derive the SZ flux-temperature relation within the tem- 
perature range of 4 keV < Tx < 11 keV. and compare it with other numerical/observational 
studies. While our findings are consistent with other SZ observations, the range of cluster 
temperatures covered by numerical simulations is too low to permit any meaningful compar- 
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ison. 

Finally, after marginalizing over gas fraction and point source contaminations, we could 
constrain the average dark matter halo concentration parameter of clusters with T^> 5 keV 

to Cyir — 3.4 Q9' 

Turning to the future prospects, in the short term, the use of WMAP 4-year maps 
should decrease our errors on /gas and concentration by a factor of 1.5 — 2. Fig. (5) shows 
the estimated number of observed clusters with z < 0.2, within different temperature bins, 
which is compared to the sample with measured temperatures used for this analyses. This 
shows that, even within the range of temperatures and redshifts resolved by WMAP, there are 
many more clusters that can be included in the analyses, if their temperatures are measured. 
Therefore, adding more clusters with observed X-ray temperatures to our catalog, e.g., from 
current and future Chandra observations, will increase the significance of our constraints. 

In the long run, this work serves as an example to show the power of combining extra- 
galactic observations in different frequencies in putting independent statistical constraints 
on the physical models of systems under study. More specifically, we have demonstrated that 
despite their low resolution, through combination with X-ray data, WMAP all-sky maps are 
capable of constraining the ICM physics at a comparable level of accuracy to their much 
higher- resolution pointed counterparts (e.g., OVRO and BIMA; see Grego et al. 2001) ^ 
Therefore, we predict that the scientific outcome of future high resolution CMB/SZ observa- 
tions (e.g., see Bond 2002) will be greatly enhanced through direct combination/correlation 
with a wide-angle deep X-ray survey, such as NASA's proposed Dark Universe Observatory 
(DUO) ^. It is needless to say that more accurate observational constraints on the ICM 
need to be supplemented by more sophisticated theoretical models which should be achieved 
through large and high-resolution numerical simulations. Such simulations are yet to be 
realized (see §5.3). 

We are grateful to Bruce Draine, Jim Gunn, Lyman Page, Matias Zaldarriaga, and in 
particular to David Spergel for their comments on the manuscript. The results in this paper 
have been derived using the HEALPix (Gorski, Hivon, & Wandelt 1999) package. 



^However, note that the errors in Grego et al. (2001) is dominated by X-ray measurements, while the SZ 
detections are significant in each cluster. This is exactly opposite to the case analyzed in this paper. 

^ http : / /duo . gsf c . nasa. gov 
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A. An Analytic Model of the Intra-Cluster Medium 

Numerical simulations indicate that the spherically averaged density distribution of 
dark matter, which also dominates the gravitational potential of galaxy clusters, may be 
well approximated by an NFW profile (Navarro, Prenk, & White 1997, hereafter NPW) 

p{r) = . I .f^\ — ' fo^^ ^ < ^vir = Cvir r„ (Al) 
(r/r,)(l + r/rj^ 

where ps and Vg are constants, and Cyir, the so-called concentration parameter, is the ratio of 
ryir to Tg- The virial radius, Tyir, is defined as the boundary of the relaxed structure, generally 
assumed to be the radius of the sphere with an overdensity of A ~ 200 with respect to the 
critical density of the universe 



/47rr3 \ /3/72\ /-Cvirrs 



ln(l + Cvir) - 



1 + Cvi, 



(A2) 



Thus, fixing the mass of the cluster, Mvir, and the concentration parameter, Cvir, at a given 
redshift (which sets the critical density for given cosmology) fixes the dark matter profile {ps 
and Ts), and the associated gravitational potential 

r ln(l + Cvir) - Cvir/(1 + Cvir) 

To model the distribution of the diffuse gas in the Intra-Cluster Medium (ICM), fol- 
lowing Suto ct al. (1998), wc assume that the gas follows a polytropic relation, i.e. Pg{r) = 
Pg{r)T{r)/ (nrrip) oc [Pgir)]^"^ , and that it satisfies Hydrostatic equilibrium in the NPW 
potential, which reduces to 

|:[r(r)+/.m,(l-7eff>(r)]=0. (A4) 



This preprint was prepared with the A AS macros v5.2. 
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Here, Pg{r), Pg{r), and T{r) are gas pressure, density, and temperature respectively, while 7eff 
is the effective polytropic index of the gas and rUp is the proton mass, fi = 4/(3 + 5X) ~ 0.59 
is the mean molecular weight for a cosmic hydrogen abundance of X ~ 0.76 . 

In order to integrate equation (A4), we need to set the boundary condition for T{r). 
Assuming an accretion shock at Tyir (Voit et al. 2003), which causes the cold infaUing gas 

to come to stop, the gas temperature behind the shock should be T(rvir) ~ In 
the spherical collapse model (Gunn, & Gott 1972) Tvir — 0.5 x r^, the turn-around radius, 
which implies v^^ = GMyir/ryir, and respectively fixes T(r) 

^ ^ -(1 - [0(r) - 0(r..)] + (A5) 

Note that, as our SZ signal is dominated by the most massive clusters (see §5), it 
is a fair approximation to neglect the non-gravitational heating/cooling processes, which 
only become significant for smaller clusters (e.g., Voit et al. 2003). Both simulations and 
observations seem to indicate that 7eff ~ 1.2 (e.g., FRBOl, Voit et al. 2003; Borgani et al. 
2004, and refernces therein), and thus, for the rest of our analyses, we will use this value. 
Also, simulations only predict a weak mass dependence for the concentration parameter, Cyir 
(e.g., NFW; Eke, Navarro, & Steinmetz 2001; Dolag et al. 2004), and thus, for simphcity, we 
assume a mass-independent value of Cyir- In §5 we explore the sensitivity of our observation 

on Cyir- 

Combining equations (A2),(A3), and (A5), we arrive at: 

nr) = /^^^^ /(,/,,,^; e,,„ 7,ff), (A6) 



where 

ln(l -I- CyirX)/x — ln(l -I- Cyir) 



/(x; Cyir, 7eff) = ^ + (1 ~ %ff > 



ln(l + Cyir) - Cyir/(1 + Cyir) 



and 

'2 GMy ^ ^'^ 



(A7) 



(A8) 



The polytropic relation can be used to obtain the ICM gas density, Pg{r) oc [T'(r)]^/('''«ff ^\ 
which yields 

Pgi''^) = /gas gir/r^ir] Cyir, Teff), (A9) 

where 

/ \ [/(^j Cyir; Teflf)] ^^'^"^ ^ /\-ir\\ 

g{x; Cyir, 7eff) = -zf— , (AlO) 

/o[/(y;Cvir,7efr)]V(Teff-i) y^dy 
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and /gas is the fraction of total mass in the ICM gas. 

Now, let us obtain the observable quantities that are relevant to our study. Similar 
to previous works (e.g., Suto et al. 1998; Komatsu & Seljak 2001), we approximate the 
observable X-ray temperature of clusters as the X-ray emission weighted gas temperature 

jT'/\r)pl{r) r'dr _ ^GM^,,m, pl\x)g\x) xHx 

- fT^Wgir) r^dr r™ r'\x)g\x) x^dx ^ ^ 

which should be contrasted with the virial (mass-weighted) temperature 

_ jT{r)pg{r) r^dr _ i^GM^i^rUp f{x)g{x) x'^dx 

Jp,{r)r^dr ' r., J,' 9{x) x^dx ' ^ ^ 

As an example, for a nominal value of Cyir = 3, we find 
while, for Cyir = 5 



We note that these relations are consistent, at the 10% level, with the predictions of the 
universal gas profile model by Komatsu & Seljak (2001). Observations of X-ray mass- 
temperature relation are often expressed in terms of M500, i.e. the mass enclosed in- 
side the sphere with A = 500 (see equation A2). For Cyir = 3(5), Equation (A2) gives 
Msoo/Myi, = M500/M200 = 0.648(0.722) yielding 



M500 
corresponding to 



/ rp \ 3/2 

cvi, = 3 or 5) = (4.0 or 3.7) x lO'^h^^' M© ( 1 , (A15) 



/ y \ 3/2 

M2oo(cvir = 3 or 5) = (6.2 or 5.1) x 10'%o' M© ( ) , (A16) 

where /170 = h/0.7. Borgani et al. (2004) argue that, on average, the beta model polytropic 
{(3^) mass estimates overestimate the normalization of the observed M-T relation by about 
30% at A = 500. Taking this into account, we notice that the normalization of the M-T 
relation in our model is consistent with the observations for hot clusters (Tx > 2 — 3keV), at 
the ~ 10% level (e.g., see Table 3 in Arnaud et al. 2005). Therefore, we will use Equation 
(All) to relate the observed X-ray temperatures of our clusters to their virial masses and 
radii. The expected systematic error in this conversion will be at the 10% level in mass 
estimates (or 3% in virial radii). 
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Fig. 1. — The distribution of cluster redshifts and virial radii (estimated from X-ray tem- 
perature, assuming Cvir = 5; see §3). The three upper hues show the resolution of WMAP 
bands (associated with the radius of the disk with the same effective area as the detector 
beams; see Page et al. 2003), while the lower line shows the physical radius of the 1 degree 
circle at the cluster redshift. 
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Fig. 2. — The binned values of gas fractions (solid circles/errorbars), inferred from our SZ 
measurements (assuming c^ij. = 3 or Cyij. = 5; Table 2). The solid triangles are the X-ray 
estimates for /gas^ (see the text), while the region enclosed by the dashed lines is the upper 
limit from the WMAP concordance model (Spergel et al. 2003). 
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Fig. 3. — The binned values of the average microwave luminosity per galaxy, Lps, assuming 
a flat spectrum. 
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Fig. 4. — The 68% and 95% likelihood contours in the fg^sh — c^ij- plane for clusters with 
Tx > 5 keV. The solid contours show the result from our SZ analysis, as described in the 
text. The dotted contours are for the ACDM concordance model, where fgash is from Spergel 
et al. (2003), while the concentration parameter and its uncertainty are obtained through 
the formalism developed in Afshordi & Cen (2002). 
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Fig. 5. — Expected cluster abundances in the nearby {z < 0.2) universe over 4.24 str (upper 
histogram), compared to the clusters used in our analysis (lower histogram). The expected 
abundance is estimated for an X-ray survey with flux limit and sky coverage similar to the 
REFLEX survey, using the concordant cosmological parameters. The comparison suggests 
that there are many more nearby clusters that can be added to our large sample, once their 
temperature is measured. 



